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Abstract. Earlier this year Chan extended the low-density series for the hard-squares 
partition function k(z) to 92 terms. Here we analyse this extended series focusing on 
the behaviour at the dominant singularity Zd which lies on on the negative fugacity 
axis. We find that the series has a confluent singularity of order 2 at Zd with exponents 
9 = 0.83333(2) and 9' = 1.6676(3). We thus confirm that the exponent 9 has the exact 
value | as observed by Dhar. 

In PP Y. Chan extended the low-density series for the hard squares partition 
function k(z) to 92 terms. In a brief analysis of the associated magnetisation series 
M(z) — 4- In k(z) Chan found that this series has a physical singularity at z c = 
3.79635(9) with exponent a = 0.0020(17). This in complete agreement with more 
accurate numerical work that has unequivocally established that the critical behaviour 
of the hard-squares model is in the Ising universality class and hence a = 0. The 
current best estimate for z c is to our knowledge z c = 3.79625517391234(4) [2J. More 
interestingly Chan gives a very accurate estimate for the dominant singularity of k(z) 
which happens to lie on the negative real axis at Zd = —0.1193388818(6) with critical 
exponent —7 = 0.171(14). Dhar [3 J has related the hard-squares model to directed-site 
animals and analysed a 42-term hard-squares series for k(z) calculated by Baxter et al 
|4] and found that the critical exponent at Zd is 9 — |. Obviously, 7 = 1 — 9 — g. 

The results of Chan are thus are quite surprising in that Zd is obtained to 10 digit 
accuracy and 7 only to 2 digit accuracy. This immediately suggests that the critical 
behaviour at Zd is more complicated than assumed in Chan's analysis. The most obvious 
complication is that the singularity at Zd contains confluent terms. In this comment we 
analyse the series for k(z) and demonstrate that this is indeed the case and we show 
that the critical point Zd is a double root. This refined analysis then allows us to obtain 
accurate estimates for the exponents at Zd, namely, 9 = 0.83333(2) and 9' = 1.6676(3), 
which obviously confirms the observation by Dhar [3] that 9 = 5/6. 

To estimate the singularities and exponents of k(z) we (as did Chan) use the 
numerical method of differential approximants [5]. We refer the interested reader to 
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Table 1. Real zeroes of Q3 and the associated exponents obtained from a homogeneous 
third order differential approximant with polynomials of degree 21. 



Zero Exponent Zero Exponent 



3.797609243403 
3.696982277961 
-0.119338882393 
-0.119338893414 

-0.120041211934 
-0.173738685995 
-0.259388829117 



2.039756900951 
7.327149084651 
0.748746994128 
2.752278507267 
2.016623831595 
8.477340013372 
1.066025488201 



-0.302260961239 
-0.379900722124 
-0.741300568486 
-0.985745094861 
-1.980070905154 
-3.832447177296 



3.431304831739 
4.228776094690 
2.959881511581 
1.581123317769 
2.502245197257 
0.541758005446 



[5] for details, and Chapter 8 of [B] for an overview of the method. Suffice to briefly say 
that a i^'th-order differential approximant to a function F(z), for which one has derived 
a series expansion, is formed by determining the coefficients in the polynomials Qi(z) 
and P(z) of order Ni and L, respective, so that the solution F(z) to the inhomogeneous 
differential equation 

j2Q t (z)(z±YF(z) = P(z) (1) 

agrees with the series coefficients of F(z). The possible singularities of the series appear 
as the zeros Z{ of the polynomial Qk{ z ) and the associated critical exponents Aj are 
obtained from the indicial equation. Note that not all roots of Qk are actual singularities 
of the underlying series. 

In tabled] we list all the real zeros of Qs(z) and the associated exponents as obtained 
from a homogeneous third order differential approximant with polynomials of degree 21. 
The exponents were calculated assuming that all the roots are distinct and hence of order 
1. We immediately notice that if the two zeros close to Za (bold-faced in the table) are 
distinct they lie incredibly close to one another. A more likely scenario is that the root 
at Zd is of order 2. If we assume this is the case and solve the resulting indicial equation 
(using the average of the two zeros for z d ) we get the exponent estimates 0.833329270 
and 1.667679940, which immediately suggests that the leading exponent is 9 = 5/6 in 
agreement with Dhar's result and that possibly the sub-dominant exponent 9' is twice 
this. The zero at ~ —0.1200 though very close to z^ is likely distinct from z&. 

In table [2] we list estimates for Zd and the two associated critical exponents obtained 
by averaging several second or third order inhomogeneous differential approximants with 
a given degree L of the inhomogeneous polynomial. The quoted error is simply the mean 
deviation of the approximants. We conclude that 

z d = -0.119338888(4), 9 = 0.83333(2), 9' = 1.6676(3). 
Clearly 9 = | confirming Dhar's result. While 9' is tantalisingly close to twice 9 our 
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Table 2. Estimates for the dominant singular point Zd and the associated exponents 
as obtained from second and third order differential approximants. L is the degree of 
the inhomogeneous polynomial. The entry for L — are estimates from homogeneous 
approximants. 



Second order approximants 



L 


Zd 






e 


& 





-0.1193388 


199(11) 


o.. 


3333312(12) 


1.66864(62) 


2 


-0.1193388 


18907(24) 


0, 


33334442(22) 


1.667787(20) 


4 


-0.1193388 


1866(21) 


0, 


8333408(38) 


1.66760(16) 


6 


-0.1193388 


1864(25) 


0. 


3333403(81) 


1.66753(23) 


8 


-0.1193388 


1868(18) 


o.< 


3333425(35) 


1.66758(15) 


10 


-0.1193388 


18754(23) 


o.. 


33334442(22) 


1.667643(20) 



Third order approximants 



L 


Zd 






e 


& 





-0.1193388 


890(18) 


0.; 


3333289(10) 


1.66777(13) 


2 


-0.1193388. 


359(24) 


o.. 


33320(30) 


1.66745(25) 


4 


-0.1193388. 


866(24) 


o.. 


33324(21) 


1.66753(28) 


6 


-0.1193388- 


349(33) 


0.; 


33318(24) 


1.66735(38) 


8 


-0.1193388847(36) 


0.; 


33318(28) 


1.66733(40) 


10 


-0.1193388. 


837(16) 


0.83348(14) 


1.66718(19) 



current best estimate does seem to exclude this possibility. 

In conclusion we have analysed the series for k(z) using differential approximants 
and found that the dominant singularity at za appears to have order 2. When this is 
taken into account the method of differential approximants is perfectly well capable of 
yielding accurate exponent estimates. In particular we confirm that 9 = | as found by 
Dhar [3]. 
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